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ABSTRACT 

We have compiled the emission-hne fluxes of O I A8446, O I A11287, and the near-infrared (IR) 
Ca II triplet (A8579) observed in 11 quasars. These lines are considered to emerge from the same 
gas as do the Fe II lines in the low- ionized portion of the broad emission line region (BELR). 
The compiled quasars are distributed over wide ranges of redshift (0.06 < z < 1.08) and of 
luminosity (—29.8 < Mb < —22.1), thus providing a useful sample to investigate the line-emitting 
gas properties in various quasar environments. The measured line strengths and velocities, as 
functions of the quasar properties, are analyzed using photoionization model calculations. We 
found that the flux ratio between the Ca II triplet and O I A8446 is hardly dependent on the 
redshift or luminosity, indicating similar gas densities in the emission region from quasar to 
quasar. On the other hand, a scatter of the O I A11287/A8446 ratios appears to imply the 
diversity of the ionization parameter. These facts invoke a picture of the line-emitting gases in 
quasars that have similar densities and are located at regions exposed to various ionizing radiation 
fluxes. The observed O I line widths are found to be remarkably similar over more than 3 orders 
of magnitude in luminosity, which indicates a kinematically determined location of the emission 
region and is in clear contrast to the case of H I lines. We also argue about the dust presence in 
the emission region since the region is suggested to be located near the dust sublimation point 
at the outer edge of the BELR. 
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1. Introduction 

Active galactic nuclei (AGNs) are known to 
have strong emission lines of various ion species. 
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Among them, the Fe II emission lines are one 
of the most prominent features in the ultraviolet 
(UV) to optical spectrum of many AGNs. They 
have long been hoped to provide significant infor- 
mation about some aspects of the AGNs and their 
host environments, e.g., energy budget of the line 
emission region and the epoch of the first star for- 
mation in the host galaxies. The determination 
of the first star formation epoch is based on the 
standard theory that the iron enrichment in galax- 
ies is delayed compared to that of the a-elements, 
such as magnesium, due to their different ori- 
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gins; Type la supernovae for i ron and Type II su- 
pernovae for the g-elem ents ( Hamann fc Ferlandl 



19931 : lYoshii et al.l 119981 ) . The delay corresponds 
to the difference in hfe-times of the progenitors of 
the two types of supernovae, and is estimated to 
be 0.3 - 1 Gyr depending on the host g alaxy envi- 



ronments (jYoshii et al.lll996l : iMatteucci fc Recchi 
20011 ). Many observations have been devoted to 



the measurement of Fe II/Mg II line flux ratios 



last decade (e.g., Elston et al.lll994 


; iKawara et al. 


1996 


: Dietrich et al. 20021 2003lj 


Iwamuro et al. 


2002 
2003 


L l20'ol Freudline et al. l2003l 


Maiolino et al. 


). However, the observed Fe II/Mg II ratios 



show a large scatter, preventing a detection of any 
significant trend in the Fe abundance as a function 
of redshift. While a part of the scatter might be 
due to the difference in the intrinsic Fe/Mg abun- 
dance ratio, it is presumable that the diversity 
of the physical condition within the line-forming 
gas, affecting line emissivities, is the ma in cause 
(jVerner et al]|2003t [BTldwin et al.ll2004 ). In the 
same sense, a change of the observed Fe II/Mg II 
ratio as a function of redshift, if found, should 
be carefully examined to tell whether it reflects 
the abundance evolution or the systematic vari- 
ation of the line emissivity. Thus establishment 
of a method to probe the line-emitting gas and 
estimate its physical parameters such as density 
and incident-ionizing radiation flux has been much 
awaited. 

Unfortunately, the Fe+ atom is characterized by 
an enormous numbers of possible electronic transi- 
tions, yielding the "Fe II pseudocontinuum" often 
observed in AGN spectra, which makes analysis 
of the observations extremely difficult from both 
the observational an d theoretical viewpoints (e.g., 
Tsuzuki et al]|2006l hereafter T06). On the other 
hand, a promising approach is to use the emis- 
sion lines emitted by simple atoms in the same re- 
gion as the Fe II lines. The most potent lines are 
O I and Ca II, whose co-spatial emergence with 
Fe II is indicated by a resemblanc e of their pro- 
files (jRodriguez-Ardila et al.ll2002a) and by a cor- 
relation between the line strengths (jPerss on 1988). 
Note that it is a natural consequence of similar ion- 
ization potentials of the relevant ions, i.e., 16.2 eV 
for Fe II, 13.6 eV for O I, and 11.9 eV for Ca II. 

The first extensive study of the physical prop- 
erties of O I emitting gases in AGNs was presented 



bv lGrandil(|l980l ). who observed the strongest O I 
line, A8446, as well as other weaker O I lines in 
Seyfert 1 (Syl) galaxies. He found that O I A8446 
lacks the narrow component that characterizes 
other permitted lines, and concluded that the line 
is purely a BELR phenomenon. He also suggested 
that O I A8446 is produced by Ly/3 fluorescence, 
which was later confirmed by the observation of I 
Zw 1, the proto t ype na rro w-hne Seyfert 1 (NLSl), 



by Rudv et al. ( 1989f) . iRodriguez-Ardila et al 



(|2002bD compiled the UV and near-IR O I lines, 
namely, A1304, A8446, and A11287, in normal Syls 
and NLSls in order to investigate their flux ratios. 
They found that there must be an additional ex- 
citation mechanism for O I A8446 — besides Ly/3 
fluorescence — which they concluded is collisional 
excitation. As for the Ca II l ines, ex t ensive 
stu dies of Syls were presente d bv iPersso ^ (|l98ah 
and Ferland fc Persson ( 19891 ). We have observed 
seven quasars at redshifts up to ~1.0 for the pur- 
pose of obtaining the UV and near-IR O I and 
Ca II lines, thus extending the previous studies to 
include quasars at high redshifts. Photoionization 
model calculations were performed and compared 
with the observations, which led us to conclude 
that the O I and Ca II lines are formed in a gas 
with density tt-h = 10^^-^ - 10^^ ° cm^^, illumi- 
nated by the ionizing radiation corresponding to 
the ionization parameter of U — 10~'^ ° — 10~^'^ 
(|Matsuoka et al.ll2005l lioorl the latter is referred 
to as Paper I hereafter). 

Now that the general properties of the O I and 
Ca II emitting g whole, are thus being 

revealed, we should turn our eyes to the follow- 
ing question: how are these gas properties re- 
lated to the quasar characteristics, such as red- 
shift and luminosity? Here we present results of 
a compilation of the emission-line fluxes of O I 
A8446, O I A11287, and the near-IR Ca II triplet 
(A8498, A8542, and A8662) observed in 11 quasars. 
The quasars are distributed over wide ranges of 
redshift (0.06 < z < 1.08) and of luminosity 
(-29.8 < Mb < -22.1), thus providing a use- 
ful sample to track the line-emitting gas proper- 
ties in various quasar environments. The observa- 
tional and theoretical data used in this work are 
described in ^ and results and discussion appear 
in iJ3l Our conclusions are summarized in fj4l 
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Data 



2.1. Observations 



A quasar sample i s com piled from lMatsuoka et al 
(|2005l ). lRiffel et all (|2006l ) and Paper I. Their ob- 
servational characteristics I are summarize d in Ta- 
ble [H The sample in iRiffel et alJ (|2006h mainly 
provides quasars in the local universe (z < 0.1); 
they selected th e quasars from the Palom er bright 
quasar survey of lSchmidt &: Green (Il983h with the 
criteria of (1) being well-known studied sources in 
the optical/UV and X-ray regions, and (2) having 
the if-band magnitude lim ited to K < 12. On the 
other hand, the quasars in Matsuo ka et al. I (|2005l) 
and Paper I constitute a high-redshift sample with 
redshifts up to 1.0. The selection criteria used in 
those works are that (1) the UV spectrum is avail- 
able in the Hubble Spa ce Telescope Faint Objec t 
Spectrograph archive (jEvans fc Koratkaii 120041 ): 
and (2) O I A8446, O I A11287, and the near-IR 
Ca II triplet fall within the atmospheric windows. 
Then the brightest objects in the target list were 
observed, which have the iJ-band magnitudes of 
H < 15. 

Although the fluxes of the releva nt O I and 



Ca II lines have been measured in iRiffel et al 
( 20061 ). we re- measured them with the identical 



method to Paper I in order to remove any sys- 
tematic differences due to the different measure- 
ment strategies, especially for deblending the O I 
A8446 -I- the Ca II triplet feature. The Icr error 
of the spectra is assumed to be 1% of the con- 
tinuum level while their spectra apparently have 
the higher qualities. The measured fl uxes of O I 
A1128 7 is plotted versus those listed in IRiffel et al ' 



()2006l) m Figure [1] they are in good agreement 
with each other, which is well within the accu- 
racy needed in this work. We also re-measured 
th e emission-line f l uxes o f PG 1116-1-215 studied 
Matsuoka et al. I (|2005l) . since we did not cover 



the Ca II lines in that work. 

^ The B-band absolute magn i tudes of the quasars are taken 
from IVeron-Cettv fc VeronI ll2006h who assumed the cos- 
mological parameters of Hq = 50 km s~^ Mpc~^ and 
go = 0. In a flat cosmology with Hp = 71 km s~ ^ 
Mpc-l, n-M = 0.29, and Ua = 0.71 (e.g.. lRiess et al.ll2004r) . 
the absolute magnitudes would be systematically smaller 
than listed in Table [l] by 0.6 - 0.7 mag (see Fig. 2 in 

IVeron-Cettv fc VeronI l2006f) . However, such a systematic 
change does not affect our arguments below. 



The measured line strengths are listed in Ta- 
ble [2] as the photon number flux ratios of the 
Ca II triplet/O I A8446 and of O I A11287/0 I 
A8446, and as the rest-frame equivalent widths 
(EWs) of O I A8446. These forms of expression 
are particularly useful in discussing the line for- 
mation processes. Hereafter they are expressed as 
n(Ca II)/n(0 I A8446), O I n(A11287)/n(A8446), 
and EW (O I A8446), respectively. 

2.2. Model Calculations 

In order to interpret the observations, we utilize 
the model calculations presented in Paper I. They 
are briefly summarized below. 

The model calculations were performed in the 
framework of the photoionized BELR gas using 
the photoionization code Cloudy, version 06.02 



(jFerland et al.l 119981 ). Two shapes of the inci- 



dent continuum were adopted, which are "stan- 
dard AGN ionizing continu um" s reported in T06 



and in iKorista et al.l (|l997l hereafter K97) . How 



ever, we found little effect on the predicted O I 
and Ca II line ratios by changing the incident con- 
tinuum shape from one to another. The BELR 
gas was modeled to have a constant hydrogen 
density nn (cm~'^) and be exposed to the ion- 
izing continuum with a photon flux of $ (s~^ 
cm~^). In the calculations, the ionizing continuum 
flux was expressed with the ionization parameter 
U = $/(rtHc) where c is the speed of light. The 
gas column density was set to 10^^ cm~^ and was 
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10 



25 



cm 



changed in a range of A^h — 10 
while the chemical abundance was assumed to be 
solar. Three cases of the microturbulence were 
considered, whose velocity vtmh is 0, 10, and 100 
km s~^. 

The different continuum shapes and the micro- 
turbulent velocities were combined into four base- 
line models (Model 1 with T06 continuum and 
Vturh = km s~^. Model 2 with T06 continuum 
and i^turb = 10 km s~^. Model 3 with TOG con- 
tinuum and wtuib — 100 km s~^, and Model 4 
with K97 continuum and wturb = 100 km s~^). 
We searched in the gas physical parameter space 
of (riH, U) in each model for the parameter sets 
that reproduce the observed O I and Ca II line 
strengths, and found that all baseline models with 



(tih, U) around (10 cm 



10 



-3.0 



) best fit to 



the observations. On the other hand, the observed 
shape and EW of the Fe II UV emission could not 
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be reproduced unless the microtur bulent velocity 
of vturh > 100 km s~^ was assumed (^Baldwin et al 



20041 ). Thus we adopt Model 3 as the default base- 



line model below, while we get little change of the 
results in this work concerning the analysis of the 
O I and Ca II lines when other baseline models are 
adopted. The (nn, U) parameters that best fit to 
the observations in Model 3 are (tt-h, U) = (10^^ ° 
cm~'^, 10"^'^), which are used as the reference grid 
point below. 

3. Results and Discussion 

3.1. A Picture of the Dust-Pree Emission 
Region 

We plot the observed values of n{Ca. II)/n(0 I 
A8446), O I n(A11287)/n(A8446), and EW (O I 
A8446) as functions of the redshift and of the B- 
band absolute magnitude Mb in Figure [21 One 
of the remarkable results is found in the top pan- 
els; the Ca II/O I A8446 ratio is hardly depen- 
dent on redshift or luminosity over the plotted 
range, while the ratio is predicted to be very sen- 
sitive to the density of the line-emitting gas in 
the photoionization models. We show the model 
predictions on the n(Ca II)/n(0 I A8446) - O I 
n(A11287)/n(A8446) plane in Figure ^ (left) , as 
well as the observed values (see also Fig. 7 in Pa- 
per The gas density nn and the ionization pa- 
rameter U are changed around the reference grid 
point, (nn, U) = (lO^^.o ^^-3^ IQ-^-^), over 2 or- 
ders of magnitude in both parameters. It is clearly 
seen that the predicted n(Ca II)/n(0 I A8446) 
ratio increases monotonically with the increased 
gas density, and that all the observed values are 
marked with the density in the vicinity of the refer- 
ence point, log riH = 12.0. Thus the similar density 
of the line-emitting gases are strongly indicated 
for the quasars distributed over these redshift and 
luminosity ranges. 

The values of EW (O I A8446), both observed 



^ Note that Figure 7 in Paper I shows the predictions of 
Model 1, which is not adopted in this paper since the as- 
sumed microturbulent velocity, i^turb = km s"'^, could 
not reproduce the observed Fe II UV emissions (see ^2.2| I. 
However, Model 3 adopted in this work predict very sim- 
ilar results to those shown in Figure 7 regarding the O I 
and Ca II emissions, while the whole pattern of contour 
is slightly {~ 0.5 dex) shifted to the high-density regime; 
the best-fit parameters in Model 1 are (riH, U) = (10^^'^ 
cm-3, 10-^°). 



and calculated with the density of log rin = 12.0, 
are shown in Figure [3] (right) as a function of the 
O I n(A11287)/ n(A8446) ratio. It shows that the 
predictions of EWs are also consistent with the 
observed data when log nn = 12.0 and a cover- 
ing fraction (cf) of the line-emitting gas as seen 
from the central energy source of 0.2 - 0.5 are 
assumed. Note that the covering fraction could 
be much smaller if we assumed oxygen overabun- 
dance relative to the solar value. 

On the other hand, a scatter of the observed 
data in the O I n(A11287)/ n(A8446) axis seems 
to be related to the diversity of the ionization pa- 
rameter (Fig. [HI left). Note that the diversity 
of other parameters, such as microturbulent ve- 
locity, gas column density, and chemical compo- 
sition, could not explain this diagram since they 
significantly alter the n(Ca II)/n(0 I A8446) ratio, 
rather than O I n(A11287)/ n(A8446), and thus 
they are unable to explain the observed similarity 
of the former ratios (Paper I). It is also quite un- 
likely that the diversity of these parameter values 
is balanced out by the fine-tuned density in such a 
way that the n(Ca II)/n(0 I A8446) ratio is always 
kept to be ~1.0, unless these lines are the dom- 
inant heating or cooling sources of the emission 
region. As with the n(Ca II)/n(0 I A8446) ratio, 
O I n(A11287)/ n(A8446) is not clearly dependent 
on the redshift or luminosity (Fig. [2l middle left 
and middle right). 

The above arguments invoke a picture of the 
line-emitting gases in quasars that have similar 
densities and are located at regions exposed to var- 
ious ionizing radiation fluxes. It would be a conse- 
quence of the difference in distance to the central 
continuum source and/or in the intrinsic luminos- 
ity of the quasars. Note that it is in clear contrast 
to the well-studied case of 11/3, whose emission 
regions in AGNs are known to be characterized 
by similar ionization parameters. In fact, rever- 
beration mapping results for 11(3 show the emis- 
sion region size (r) - luminosity (L) relation of 
r oc L'^-^, which is consiste nt with the constant 
ionization p arameter regime ([Peterson et al.ll2002t 



Bentz et al.ll2006i) . Such a situation has long been 



expected in order to account for the remarkably 
similar AGN spectra over a broad range of lumi- 
nosity, and was incorporated into the locally op- 
timally emitting clou d (LOG) model suggested by 
Baldwin et"aD (|l995f ). that is, that the BELR is 
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composed of gas with widely distributed physi- 
cal parameters and each emission line arises from 
its preferable environment. On the other hand, 
the case in O I and Ca II lines apparently indi- 
cates that the location of the emission region is 
not radiation-selected. 

In line with the above arguments, we found a 
clear difference of the velocity-luminosity relation 
between O I and H/3; the measured O I line widths 
are plotted versus Mb in Figure [H which shows 
that the O I line widths are remarkably similar, 
concentrated around 1500 - 2000 km s~^, over 
more than 3 orders of magnitude in the i?-band 
luminosity. On the other hand, those for the hy- 
drogen Balmer lines usually have a large scatter as 
shown by, e.g.. lKaspi et al. ( 2000f) : their sample of 
34 AGNs, spanning over 4 orders of magnitude in 
continuum luminosity, has the line widths of 1000 
- 10,000 km s~ ^. Such a trend is also indicated by 
Persson (Il988h . who reported that while the cor- 



relation between FWHM (O I) and FWHM (H/3) 
is good for the small FWHM regime, the O I lines 
grow systematica lly narrower than H/3 at large line 
width. iRodriguez-Ardila et al. ( 2002a ) conducted 
a detailed study of the near-IR emission line pro- 
files in NLSls, and found that O I, Ca II, and Fe II 
lines are systematically narrower than the broad 
components of other low-ionization lines such as 
hydrogen Paschen lines and He I A10830. They 
also argued that these lines are produced in the 
outermost portion of the BELR, since their widths 
are just slightly broader than those of [S III] A9531 
which they assumed is formed in the inner portion 
of the narrow emission line region (NELR). While 
the scattered widths of H I lines could be inter- 
preted as a consequence of the radiation-selected 
locations of the emitting gases, regardless of the 
gas kinematics, the remarkable similarity of the 
O I line widths might imply the kinematically de- 
termined emission regions. In such a situation, the 
diversity of the ionization parameters, as discussed 
above, would b e an inevitable result. As stated by 
IPerssonI (1988), there is clearly interesting infor- 
mation which could be deduced from studies of 
the O I and H I line profiles; especially, reverbera- 
tion mapping of these O I and Ca II lines would be 
a powerful tool to reveal the underlying physics. 



3.2. A Picture of the Dusty Emission Re- 
gion 

The widely-accepted theory of the BELR de- 
scribes its outer edge, where the O I and Ca II 
emission lines are likely t o be formed, set b'v 



emission lines are likely t o be tormed, set by 
the dust sublimatio n (e.g., Laor fc Draind 19931 
Netzer fc Laoi1l993h . If we accept this picture, the 



dust grains are possibly mixed in the line-emitting 
gas and suppress the Ca II emission through the 
substantial Ca depletion. Such a situation is in 
fact reported for the NELR by the absence or sig- 
nificant weakness of the obse rved [Ca II] A7291 
line (Kingdon et al. 1995; Vill ar-Martfn fc Binette 
ri997t ). iFerguson et all (|l997t) presented the LOC 
model calculations of the narrow emission lines 
and argued that Ca is depleted relative to the solar 
value by factors of 3 - 160. 

It is hard to see an evidence of the dust pres- 
ence in the emission region from our results, since 
Figure [3] appears to show that our dust-free mod- 
els successfully reproduce the observations. How- 
ever, it is noteworthy that it is quite difficult to ac- 
count for the observed data at O I n{\ 11287)/n(A 
8446) < 0.4 by the models with log nn = 12.0. 
The problem is that such small values of O I n{\ 
11287)/n(A 8446) could only be reproduced with 
the higher densities than log n-n — 12.0 so that 
O I A8446 emission is exclusively enhanced by the 
collisional excitation, while such a dense gas pro- 
duces intense Ca II emission that is much stronger 
than observed. One can clearly see this trend in 
Figure [3] (Ze/i(). If we assumed significant Ca deple- 
tion in the line-emitting gas, these difficulties are 
naturally resolved since it significantly suppresses 
the otherwise intense Ca II emissions. For exam- 
ple, the data point representing (log nu, log U) 
— (13.0, —2.5) in Figure [3] {left) could provide a 
plausible model for the observations at O I n{\ 
11287)/n(A 8446) < 0.4 if the Ca II emission is 
suppressed by a factor of a few times loH 

The dust grains present in the line-emitting 
gas might also give the natural explanation to 
the observed lack of some UV emission lines rel- 
ative to their optical or near-IR counterparts. 



^ The predicted EW for the model with (log riHi log U) = 
(13.0, -2.5) is EW (O I A8446) = 30 A, which is consistent 
with the observations if rather large covering factors of 0.5 
— 1.0 and/or oxygen overabundance relative to the solar 
value are assumed. 
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Ferland fc Persson ( 19891 ) mentioned the possibil- 
ity of the dust survival in the BELR gas in or- 
der to explain the extreme weakness of the ob- 
served Ca II A3934 and A3639 lines relative to 
the near-IR triplet. At least a part of the long- 
standing Fe II UV/opt problem, in which the ob- 
served ratios of Fe II UV flux to the optical flux 
fall far below the photoionization model predic- 
tions (see, e.g.. [Baldwin et al.ll2004[ ). could also be 
explained. However, it should be noted that the 
dust would affect the line formation processes in 
a very complicated manner, which should be pre- 
cisely addressed when discussing the specific lines; 
for the O I and Ca II lines, one of the most appar- 
ent effects as well as the Ca depletion would be 
the destruction of Ly/3 photons which otherwise 
excite O I atoms, thus the suppression of the O I 
emissions. The resultant line flux ratios could be 
much different from those derived from the simple 
speculations. 

4. Summary 

We have compiled the emission-line fluxes of O I 
A8446, O I A11287, and the near-IR Ca II triplet 
(A8579) observed in 11 quasars. The quasars 
are distributed over wide ranges of redshift (0.06 
< z < 1.08) and of luminosity (-29.8 < Mb < 
—22.1), thus providing a useful sample to track 
the line-emitting gas properties in various quasar 
environments. The measured line strengths and 
velocities, as functions of the quasar properties, 
were analyzed with photoionization model calcu- 
lations. Our findings and conclusions are as fol- 
lows: 

1. There is no sign of a significant change in the 
flux ratios of the Ca II triplet and O I A8446 over 
the redshift and luminosity ranges studied here. 
It strongly indicates similar gas densities in the 
line-emission region from quasar to quasar. 

2. The observed scatter of the O I A11287/A8446 
ratios appears to be related to the diversity of the 
ionization parameter, while the ratio is not clearly 
dependent on the redshift or luminosity. Com- 
bined with the similarity of the Ca II/O I A8446 
ratios, it invokes the picture of the line-emitting 
gases in quasars that have similar densities and 
are located at regions exposed to various ionizing 
radiation fluxes. 

3. The O I line widths are remarkably similar 



from quasar to quasar over more than 3 orders of 
magnitude in luminosity. It might imply a kine- 
matically determined location of the line emission 
region and is in clear contrast to the case of H I 
lines, whose emission region is considered to be 
radiation-selected. 

4. If we accept that the O I and Ca II emission 
lines are formed at the outer edge of the BELR and 
that the outer edge is set by the dust sublimation, 
the line-emitting gas is possibly mixed with the 
dust grains. In fact such a situation may better 
reproduce the observations than the dust-free case 
through the significant Ca depiction. 
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Table 1: Sample Characteristics 



Quasar 


Redshift 




Reference'' 


QSO B0850+440 


0.514 


-25.2 


3 


PG 1116+215 


0.176 


-25.3 


2 


PG 1126-041 


0.060 


-22.8 


1 


QSO J1139-1350 


0.560 


-26.3 


3 


PG 1148+549 


0.978 


-28.0 


3 


3C 273 


0.158 


-26.9 


3 


PG 1415+451 


0.114 


-22.1 


1 


PG 1448+273 


0.065 


-23.0 


1 


PG 1519+226 


0.137 


-22.9 


1 


PG 1612+261 


0.131 


-22.6 


1 


PG 1718+481 


1.084 


-29.8 


3 



"The S-band absolute magnitude taken from Iviron- Getty fc Veror] ||2006D . 
''References — (D iRiffel et"ai] l l2006t) : (2) iMatsuoka et al.l ||2005|) : (3) Paper I 



Table 2: Observed Photon Number Flux Ratios 
and EWs of the O I and Ca II Emission Lines 



Quasar n(Ca II)/n(0 I A8446) O I n(A11287)/n(A8446) EW (O I A8446) [A] 



QSO B0850+440 


1.63 


± 


0.53 


0.62 


± 


0.15 


24.9 ± 5.9 


PG 1116+215 


1.22 


± 


0.43 


0.67 


± 


0.20 


19.8 ± 5.2 


PG 1126-041 


1.12 


± 


0.06 


0.52 


± 


0.04 


38.2 ± 1.4 


QSO J1139-1350 


0.49 


± 


0.16 


0.24 


± 


0.05 


16.0 ± 2.9 


PG 1148+549 


1.18 


± 


0.05 


0.33 


± 


0.02 


26.0 ± 0.9 


3C 273 


1.46 


± 


0.21 


0.33 


± 


0.04 


17.2 ± 1.9 


PG 1415+451 


0.87 


± 


0.12 


0.35 


± 


0.05 


23.3 ± 2.2 


PG 1448+273 


0.34 


± 


0.12 


0.54 


± 


0.08 


22.1 ± 2.5 


PG 1519+226 


0.97 


± 


0.07 


0.41 


± 


0.03 


31.9 ± 1.6 


PG 1612+261 


0.47 


± 


0.07 


0.44 


± 


0.05 


35.5 ± 2.3 


PG 1718+481 


1.11 


± 


0.22 


0.47 


± 


0.12 


7.6 ± 0.8 



Note. — The representative wavelength of 8579 A is used to convert the Ca II triplet fluxes to the photon number fluxes. 
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10 15 20 

01 Flux (Riffel + 06) 

Fig. 1.— The fluxes of O I A11287 measured 
in this work are plotted versus those Usted in 
Riffel et alJ (|2006l) . The units of both axes are 
10""'^^ ergs cm~^ s~^. The fluxes measured in PG 
1126—041 (at the top-right corner) was multiphed 
by 0.5 in both axes in order to improve the visibil- 
ity. A dotted line represents the loca tions where 
our O I fluxes are identical to those bv lRiffel et al 
()2006l) . 
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Fig. 2. — The observed values of n(Ca II)/n(0 I 
A8446), O I n(A11287)/n(A8446), and EW (O I 
A8446) are plotted as functions of the redshift 
(left) and of the i3-band absolute magnitude 
(right). 
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Fig. 3. — Observed and theoretical values of 
n(Ca II)/n(0 I A8446) (left) and EW (O I A8446) 
(right) versus O I n(A11287)/n(A8446). Filled 
circles represent the observations while diamonds 
represent the model predictions. The models with 
the same gas density, log tt-h — 11.0, 11.5, 12.0, 
12.5, and 13.0 (left), or with the same covering 
factor, cf = 0.2, 0.5, and 1.0 with log nn = 12.0 
(right), are connected with dashed hues. Gray 
scale filling the diamonds distinguishes different 
values of the ionization parameter U as indicated 
at the bottom-left corner of the left panel. 
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Fig. 4.— The observed line widths of O I A11287 
versus the B-band absolute magnitudes. 



9 



